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ABSTRACT 



We construct a magnetic helicity conserving dynamo theory which incorpo- 
rates a calculated magnetic helicity current. In this model the fluid helicity plays 
a small role in large scale magnetic field generation. Instead, the dynamo process 
is dominated by a new quantity, derived from asymmetries in the second deriva- 
tive of the velocity correlation function, closely related to the 'twist and fold' 
dynamo model. The turbulent damping term is, as expected, almost unchanged. 
Numerical simulations with a spatially constant fluid helicity and vanishing re- 
sistivity are not expected to generate large scale fields in equipartition with the 
turbulent energy density. The prospects for driving a fast dynamo under these 
circumstances are uncertain, but if it is possible, then the field must be largely 
force-free. On the other hand, there is an efficient analog to the a — Q dynamo. 
Systems whose turbulence is driven by some anisotropic local instability in a 
shearing flow, like real stars and accretion disks, and some computer simula- 
tions, may successfully drive the generation of strong large scale magnetic fields, 
provided that d r fl{dgv z U0) > 0. We show that this criterion is usually satisfied. 
Such dynamos will include a persistent, spatially coherent vertical magnetic he- 
licity current with the same sign as — d r Q, that is, positive for an accretion disk 
and negative for the Sun. We comment on the role of random magnetic helicity 
currents in storing turbulent energy in a disordered magnetic field, which will 
generate an equipartition, disordered field in a turbulent medium, and also a 
declining long wavelength tail to the power spectrum. As a result, calculations 
of the galactic 'seed' field are largely irrelevant. 
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1. Introduction 

Astrophysical dynamos are usually discussed in terms of mean-field dynamo theory 
(see Moffatt 1978; Parker 1979; Krause & Radler 1980). This typically involves several 
different assumptions. First, one assumes a dynamical separation between large scale 
and small scale (i.e. turbulent eddy scale) fields. Second, the latter are assumed to 
be characterized by the turbulent velocity field acting on the large scale field without 
systematic velocity field effects due to the large scale field (no 'back-reaction'). Third, 
the large scale field is assumed to be smoothed by turbulent diffusion. (This is follows 
from the preceding point and the assumption that reconnection is fast and efficient.) 
Finally, the velocity field properties are typically prescribed a priori, even when they are 
driven by magnetic field instabilities. In this theory, the growth rate of the large scale 
magnetic field is driven by the fluid helicity. Only the first point is a necessary part of 
dynamo theory. In fact, all other parts of mean-field dynamo theory have serious problems 
(see, for example, Cattaneo & Vainshtein 1991; Parker 1992; Kulsrud & Anderson 1992; 
Gruzinov & Diamond 1994; Gruzinov & Diamond 1996; Cattaneo & Hughes 1996; 
Brandenburg 2000). Obviously, the evidence for widespread dynamo activity in stars 
suggests that the problems here lie in our understanding of physics rather than suggesting 
that fast dynamos are impossible. 

We will not address the issue of fast reconnection here. There is a large body of evidence 
indicating that magnetic fields in astrophysical plasmas can reconnect at speeds approaching 
the Alfven speed (see Dere 1996; Innes et al. 1997, and references contained therein). This 
may be understood in terms of field line stochasticity (Lazarian & Vishniac 1999) although 
there are competing explanations (Petschek 1957; Shay et al. 1998; see also Kulsrud 2000 
for a modified version of Petschek's theory). In addition, the tendency to prescribe the 
small scale velocity field can be seen as a largely formal problem. Even when the turbulence 
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is driven by magnetic instabilities, e.g. the Balbus-Hawley instability in accretion disks 
(Balbus & Hawley 1991; Hawley & Balbus 1991), as long as the eddy scale is much smaller 
than the large scale magnetic field scale the problem can still be divided into two conceptual 
steps: solving for the properties of the small scale turbulence, and understanding the 
generation of the large scale field. 

Here we will focus on a particular aspect of the problem of magnetic back-reaction. 
We begin by noting that magnetic helicity, A • B, is strictly conserved for resistive MHD, 
as we take the resistivity to zero. In contrast, while fluid helicity v • ou, is also conserved 
in the limit of zero viscosity, the dissipation of fluid helicity at small scales does not 
vanish as v — > 0, so that fluid helicity is not even approximately conserved in any real 
turbulent system. Unfortunately, standard mean-field dynamo theory does not reflect this 
conservation law. There must be a correlated backreaction that enforces it, and this effect 
is left out of the standard theory. 

Numerical simulations seem to confirm the notion that there is a serious problem 
with mean field dynamo theory. Computer simulations of dynamos can be divided 
into two classes. There are simulations in which some local instability (convection, the 
Balbus-Hawley instability etc.) is allowed to operate, and there are simulations in which 
the turbulence is driven externally, usually in such a way as to guarantee the presence of 
a net fluid helicity. The former simulations are often successful at generating large scale 
magnetic fields whose energy density is at least as great as the turbulent energy density 
(e.g. Hawley & Balbus 1992; Glatzmaier & Roberts 1995; Brandenburg et al. 1995). The 
latter are less successful, in the sense that the energy density of the large scale magnetic 
field is often relatively modest (e.g. Meneguzzi et al. 1981; Balsara 2000). In recent years 
numerical simulations (Cattaneo & Hughes 1996; Hughes et al. 1996; Brandenburg 2000) 
have been performed to test the validity of mean field dynamo theory and the role of 
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magnetic helicity conservation in suppressing dynamo effects. These calculations have used 
a closed box containing with some sort of forced heliacal turbulence. Cattaneo & Hughes 
(1996) were able to show a strong suppression of the turbulent dynamo in the limit of small 
resistivity. Hughes, Cattaneo and Kim (1996) did find dynamo action, defined as the growth 
of the total magnetic energy density, but it did not depend on the fluid helicity. Their eddy 
scale was very close to their box size, so a clean separation between eddy scale and large 
magnetic fields was not possible. Finally, Brandenburg (2000) produced an example of a 
dynamo in a computational box, with forced heliacal turbulence. By varying the resistivity 
he was able to demonstrate a steep inverse correlation between the dynamo growth rate and 
the conductivity. Naively extrapolating to astrophysical regimes suggests that magnetic 
dynamos driven by fluid helicity would take enormous amounts of time to grow. 

This result has been anticipated by a series of analytic arguments and computational 
studies (e.g. Vainshtein and Cattaneo 1992; Gruzinov and Diamond 1994; 1996) which 
point to a suppression of the electromotive force associated with the fluid helicity. There 
have been attempts to calculate this back reaction under various approximation schemes 
(e.g. Field, Blackman and Chou 1999), but they involve an expansion in parameters that 
are typically of order unity 1 . This is not surprising, since the problem is similar to, but more 
complicated than, attempts to derive the detailed properties of fluid turbulence analytically. 
Vainshtein and Cattaneo (1992) (see also Gruzinov and Diamond 1994) have suggested, 
from fairly basic considerations, that this back reaction should be overwhelmingly strong 
as soon as the magnetic field reaches equipartition with the surrounding turbulence on the 
dissipation scale and that it should suppress the critical components of the fluid helicity 
tensor (i.e. those that contribute to an electromotive force parallel to B) on large scales. 

x The expansion of Field et al. involves a parameter which is almost exactly one in 
numerical simulations (see Cho and Vishniac 2000) 
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Since this criterion is satisfied when the large scale field is negligibly small, this looks like 
a fairly powerful anti-dynamo argument. It is consistent with the numerical work showing 
dynamo suppression. 

Here we explore the possibility that a new version of mean field dynamo theory, 
modified to explicitly incorporate magnetic helicity conservation, can be used to predict 
the evolution of large scale magnetic fields in highly conducting fluids. In section 2 of this 
paper we derive a new set of dynamo equations and apply them to simple dynamo models. 
We find that not all simple dynamos are eliminated, or even suffer reduced growth rates. In 
section 3 we discuss the role of random velocities in building a disordered magnetic field, 
and apply our results to the early evolution of the galactic magnetic field. In the final 
section of this paper we summarize our results and discuss some of the broader implications 
of this work. 

2. Magnetic Helicity Conserving Magnetic Field Evolution 

We start with the usual expression for magnetic field evolution in ideal MHD. 



where the subscript 'T' denotes the total field. In what follows, lower case letters will stand 
for fluctuation quantities so that B T = B + b. Defining the vector potential, A, in the 
usual way we can show that 




(1) 



d t (A T ■ B T ) = V • ((v x B T + V$ T ) x A T ) , 



(2) 



where $ T is defined through 



V 2 $ T = V • (v x B T ), 



(3) 



and we have used the gauge condition V • A T = 0. 
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Equation (2) guarantees global conservation of magnetic helicity. It is equally obvious 
that the gauge dependence inherent in the definition of magnetic helicity means that this 
conservation law is gauge dependent. We have chosen V • A T = 0, but any gauge condition 
that does not include eddy scale terms will give equivalent results. It is tempting to use 

d t A T = v x B T (4) 

as our gauge condition instead, so that we can drop V$t from equation (2), but this 
introduces a steady growth of magnetic helicity on eddy scales and we would no longer 
expect large scale helicity to be separately conserved. 

Mean field dynamo theory is an attempt to follow the magnetic field averaged over 
many eddies, without a detailed calculation of individual eddy scale features. To this end 
we rewrite equation (1) as 

d t B = V x £ mf , (5) 

where 

S mf = (v x b). (6) 

The usual approach is to calculate £ m f using general assumptions about the structure of the 
underlying turbulence. Assuming that velocities are correlated only when they are in the 
same direction, and that the energy is distributed roughly isotropically among the available 
modes we find that 

£ m f,i = a ijBj — €ijkDdjB k , (7) 

where 

D = ^ < v 2 > r c , (8) 

OLij = Ulm(vidjV m )T c , (9) 

and r c is a velocity correlation time. The tensor is the fluid helicity tensor. (Sometimes 
the r c is omitted from the definition.) Its trace, without the factor of r c , is the fluid helicity. 
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The difficulty with this approach is that, in general, it does not conserve magnetic 
helicity. If we rederive equation (2) following only the dynamics of the large scale fields, as 
defined in equation (5), we find that 

d t (A-B) = A • V x S mf + B • (£ mf - V$) (10) 
= 2B • £ m f + V • ((S m f + V$) x A) . 

On the other hand, we can calculate the evolution of the average magnetic helicity density 
by averaging equation (2), over eddy scales. We obtain 

d t A ■ B = V • ((S mf + V$) x A + ((e m/ + V0) x a)) . (11) 

Here we have assumed that 

(A T • B T ) = A • B, (12) 

which will be true for an eddy scale smaller than the typical field scale by a factor of at 
least b 2 /B 2 and assuming efficient transfer of magnetic helicity between scales. The latter 
assumption is examined at length in the appendix. We will consider the effects of a non-zero 
eddy scale in the third section of this paper. 

The second contribution to the magnetic helicity current is produced by correlations 
on eddy scales. We will refer to it as the anomalous magnetic helicity current, J# (with 
an implied minus sign, since it appears on the right hand side of equation (11)). Any 
self-consistent evolution equation for the large scale magnetic field has to satisfy both of 
the preceding equations. This does not uniquely specify the correct expression for S m f, but 
if we further specify that this expression should depend only on the local value of B, the 
helicity current source term in equation (11), and the conventional expression for S m f then 
we are forced to choose 

£mf = £mf,± ~ 7T^ V ' J #> ( 13 ) 
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where £ m /,j_ is the electromotive force perpendicular to the direction of the large scale 
magnetic field. We will assume that this quantity can be calculated in the conventional 
manner. In any case, the critical element for a successful mean field dynamo is the 
component of the electromotive force along the direction of the large scale magnetic field, 
and this is uniquely specified. The form of £ m f ■ B given in equation (13) was first proposed 
by Bhattacharjee and Hameiri (1986), who further constrained the form of Jh by requiring 
that the energy dissipated by small scale currents be balanced locally by the energy put into 
the large scale field. We will not invoke this constraint here, since in typical astrophysical 
systems the energy flow through the large scale magnetic field is a small fraction of the 
energy flow through the turbulent cascade, and does not need to be balanced locally with 
ohmic dissipation to conserve energy. While this work has since been extended to explore 
the weak and strong magnetic field limits, and the role of hyperresistivity in magnetic 
field evolution (Bhattacharjee and Yuan 1995), it has not yet been fully applied to the 
generation of astrophysical magnetic fields. We note that equation (13) also has the effect of 
eliminating the turbulent dissipation of currents aligned with the large scale magnetic field. 
The implication is that force-free magnetic fields are not subject to turbulent dissipation, 
but all others are. Gruzinov and Diamond (1994) concluded from this that the lifetimes of 
large scale astrophysical magnetic fields against dissipation are not qualitatively altered by 
magnetic helicity conservation. Here we note that this argument leaves open the possibility 
that force-free magnetic fields might be generated even by a very slow dynamo process. 
Despite this loophole, we can conclude that the fast generation of large scale magnetic fields 
depends on the generation of large scale helicity currents from eddy-scale processes. 

In order to evaluate the helicity current term, we take 

a = (e m/ - V0)r c , (14) 
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where, as before, r c is the eddy correlation time. Then equation (13) becomes 

4/ = ^/,i + ^V-((V^xe m/ )r c ), (15) 

(This equation can also be recovered from equation (20) in Bhattacharjee and Yuan (1995) 
in the limit where the resistivity goes to zero and the eddy size is assumed to be much 
smaller than the scale of the mean magnetic field.) The source term in equation (15) can 
be evaluated explicitly by writing <fi in terms of the Fourier transform of e m /. We find 

(V0xe m/ ) ^(/(j(I-tt)«.e m/ (x + r))lrxe m/ (x)) (16) 

= / fee ijfc 9 fc ^(e m/j -(x)e m/ ,,(x + r)), 

where the partial derivatives are with respect to the components of r. This expression can 
be further simplified by using the first-order smoothing approximation 

e mf = v x (B + b) - (v x b) « v x B (17) 

and ignoring the gradient of the mean magnetic field over eddy scales 2 . The use of the 
first-order smoothing approximation is somewhat controversial, but in this context it does 
not violate any basic conservation laws. This expression neglects the random component of 
e m f. Its effects are discussed in the next section. We find that J H is 

Jh = -(V0 x e m/ )r c = - / — ei nm B k Bid k d m {vi{^)v n {-x. + r))r c (18) 

We can rephrase this quantity in terms of correlations between derivatives of the 
velocity field if we are willing to accept a certain loss of precision. Expanding t> n ((x + r) 
around r = we find that equation (18) becomes 

Jh ~ -A 2 e ZnmJ B fcJ B;(v i (x)a m a fc v n (x))r c , (19) 



2 Turbulent dissipation effects come from including the field gradient terms when their 
coefficients are non-zero for isotropic turbulence. No such terms appear here. 
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where A c is some sort of angle-averaged eddy correlation length, roughly the geometric mean 
of the two largest perpendicular correlation lengths. If we integrate this by parts we obtain 



Clearly there is no simple relation between this and the fluid helicity. In particular, a 
computer simulation which is tailored to give a uniform non-zero fluid helicity will not 
normally produce a non-zero magnetic helicity current. However, there is an attractive 
physical interpretation to equation (20). The anomalous magnetic helicity current is 
proportional to the correlated product of the gradient of the velocity field along the 
magnetic field lines (a 'fold') and the vorticity along the magnetic field lines (a rotation). 
If one repeats this to form a complete dynamo cycle, it is obviously closely related to the 
'twist and fold' dynamo model first proposed by Vainshtein and Zel'dovich (1972). 

As a simple example we can consider the generation of a magnetic field in a differentially 
rotating flow, Q(r) oc r~ q . Galaxies, accretion disks, and stars are all instances of this case 
(although the use of cylindrical geometry is a bit suspect for stars). This is conventionally 
explained as an example of an L a — dynamo, where 



^<(B.c)(B-V)v)>r c . 



(20) 



d t B r = —d z 8 m ffi, 



(21) 



and 



d t B e = -qttB r + d z (Dd z B e ) . 



(22) 



For B r a small, fixed fraction of Be we can write equation (21) as 




(23) 



For homogeneous anisotropic turbulence this becomes 



[d t - Dd 2 z ] 2 B r = 4qQX 2 c r c (uj e -deV z )d 2 z B r . 



(24) 
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In other words, a successful dynamo in this system requires that qd$v z be negatively 
correlated with uig and that the growth rate exceed the turbulent diffusion rate. 

It is useful to compare this result with the result of a conventional, fluid helicity driven, 
l a — JT dynamo. In this case, ignoring dissipation, we get a growth rate 

r~(f!^) 1/2 , (25) 

where L is a large vertical scale associated with the structure of the magnetic field and/or 
the structure of the disk or star. Since agg is odd for a reversal of any coordinate direction, 
we need symmetry breaking in all three directions before we expect a non-zero value. 
The differential shear, and the resulting Coriolis forces, breaks symmetry in the r and 9 
directions, while ensuring symmetry under the transformation (f, 9) — > (— f, —9) (which also 
leaves agg unchanged) . We need to appeal to some kind of background vertical structure to 
provide symmetry breaking in the z direction before we can expect a nonzero agg. If we 
assume that the typical eddy velocity is ~ A c /r c , then agg ~ v \ c (Qr c )/L ~ \ 2 c Vl/L and 

r ~ oh, (26) 

for the weak shearing limit, Qt c < 1. In contrast (ujg^dgv z ) is unchanged for z — > —z, and 
will generally be non-zero for turbulence in a shearing flow. No vertical structure, aside 
from the vertical scale of the magnetic field, is necessary. Using equation (24) to estimate 
the dynamo growth rate we find that 

T ~ tt^. (27) 
Lb 

Here Lb is the vertical scale of the magnetic field, which may be considerably less than 
the background structure scale, depending on the efficiency of turbulent damping. The 
implication is that the growth rate of stellar and disk dynamos is at least as high in this 
model as it is in models which ignore magnetic helicity conservation, provided that the local 
turbulence supplies the correct sign of (ug^dgv z ). 
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How likely is this? Invoking incompressibility and integrating by parts, we can rewrite 
this quantity as: 

(uj e ^d e v z ) = (d z v z \d e v r ) - (d r v z \d e v z ) (28) 

= -(d r v z \d e v z ) - {d r v r \d v r ) - {\d e v e \d e v r ) . 

If we replace these derivatives with wavenumbers, then we have 

(oj0-d e v z ) ~ -k r k e (v 2 + v 2 z ) - k 2 (v r v e ). (29) 
r 

The effect of shear is to make k r ke positive for a positive q. In other words, the first term 
always has the correct sign for driving a dynamo. The second term may have either sign, 
but will only be significant in the presence of a strong angular momentum flux. For the 
Balbus-Hawley instability in an accretion disks, q = 3/2 and (v r ve) > 0. Consequently, both 
terms in equation (29) are negative and a dynamo is guaranteed for sufficiently large vertical 
domain size. In stars with active convection zones, like the Sun, the result is less obvious, 
but will still favor dynamo activity as long as the angular momentum distribution in the 
star is stationary, or nearly so. One striking counter-example is provided by simulations of 
the magnetic Kelvin- Helmholtz instability (Hawley, Gammie & Balbus 1996), in which they 
repeated their simulations of a zero-flux azimuthal magnetic field embedded in an accretion 
disk, but turned off the centrifugal force term. In this case the absence of a centrifugal force 
term leads to a large negative (v r vg), through turbulent mixing, and no significant dynamo 
activity was observed. 

We can understand the dynamics of this modified version of the L a — f2' dynamo by 
comparison with the standard model. In the standard picture the fluid helicity induces a 
systematic spiral twisting of the azimuthal field lines. Turbulent smoothing gives a radial 
field provided that the twisting process gives a vertical gradient in the product of the 
magnetic field strength and the pitch angle of the spiral twisting. By contrast, the process 
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considered here works through the action of the correlation (v z deue) (which is related to 
the left hand side of equation (29) through an integration by parts). This corresponds to 
twisting a bundle of field lines in one direction when the bundle is displaced in the positive 
z direction, and twisting it in the opposite sense when the z displacement is reversed. 
Turbulent smoothing gives us a net local excess of magnetic helicity when this process 
is not uniform, that is, when Jh has a non-zero divergence. Subsequently the turbulent 
smoothing of bundles of spiralling magnetic field lines proceeds as in the standard model. 
Near a maximum in \Bg\ the effect of reconnection and smoothing is to produce a coherent 
B r which, for the appropriate sign of the correlation, will reinforce the azimuthal field 
component via shearing effects. This may seem like a rather awkward substitute for the 
analogous effect in the standard model, in which azimuthal field lines are distorted into 
spirals with a coherent helicity. However, this new model involves only motions which are 
clearly realizable through a continuous distortions of the fields lines, except for the final step 
of reconnection, which involves reattachment of adjacent field lines of differing orientation. 
The standard picture involves a systematic rearrangement of the topology of the field lines 
within the bundle, where the attachments of field lines are not switched, but rather are 
allowed to slip over the surface of a plane extending through the bundle. The difference can 
be readily appreciated by anyone who has ever attempted to twist a rubber torus. This is 
the reason why the standard model involves a violation of magnetic helicity conservation, 
whereas the model explored in this paper does not. 

This particular example of a successful dynamo has a curious feature. The total 
magnetic helicity current is 



where we have dropped V$ = d z $ and £ m /,j_ since they are both zero in this model. The 




(30) 
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dominant component of the A is 

A _J Z Be , 

where k z is the vertical wavenumber of the magnetic field. Using equation (20) we can 
rewrite equation (30) as 



This is a spatially constant, vertical magnetic helicity current. Its divergence is zero, since 
the magnetic field has zero magnetic helicity, but its presence is a necessary part of the 
dynamo process. We can add oscillating pieces by making different gauge choices, but the 
large scale average is gauge invariant. In a real system this flux will appear at the vertical 
boundaries of the system. Its sign is a function of q. For an accretion disk J to t ■ z > 0. 
For a successful dynamo with q < (like the Sun), we expect a persistent negative 
magnetic helicity current. Despite the fact that the solar dynamo is undoubtedly more 
complicated than this simple example, this effect is observed (Hale 1927; Richardson 1941; 
Seehafer 1990; Webb 1992; Pevtsov et al. 1995; Rust & Kumar 1996; Canfield et al. 1999). 

Finally, while this simple model has no net accumulation of magnetic helicity, this is 
partly due to assuming perfect coupling between the long wavelength magnetic helicity 
components due to small scale structures, which we refer to here, and in the appendix, as 
h, and the long wavelength magnetic helicity components due to the large scale field itself, 
denoted H. We can estimate the size of the fluctuations in h from equation (A9), which 
implies that the fluctuations will be independent of B, and reduced from the maximum 
magnetic helicity which can be stored on the eddy scale by a factor ~ {k z \)~ l . 

On the other hand, if we consider turbulence without a large scale shear, then the 
role formerly played by the fluid helicity is played by — B~ 2 V ■ J#. While this will be 
non-zero, even for homogeneous turbulence, it is hard to imagine a situation where it can 
dominate over the turbulent dissipation term. In general, the two terms will be of the same 




(32) 
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order. The former can be made larger only at the expense of having a short wavelength 
associated with the global properties of the turbulence, producing a rapidly oscillating 
electromotive force. Moreover, the dissipation term depends on the diagonal terms in the 
Reynolds tensor, which will be at least as large as the asymmetric terms which drive J#. 
While this seems to be a strong argument against the existence of a viable analog to the 
'a 2 ' dynamo of conventional mean-field dynamo theory, there is a possible loophole. Since 
the magnetic helicity conservation constraint suppresses turbulent dissipation for force-free 
magnetic fields, it follows that a dynamo process which preferentially generates such a field 
may be possible. We will defer any further discussion of this point to a subsequent paper. 

3. Random Effects in MHD Turbulence 

We have assumed in this discussion that the individual eddies are negligibly small and 
that the only effects worth considering are those that are coherent over many eddies. In 
typical astrophysical systems, this is an exaggeration, and it is possible to generate some 
large scale field even in the absence of an effective large scale fast dynamo of the kind 
discussed above. We note that simulations of MHD turbulence usually show a magnetic 
field whose energy density is 10% or more of the turbulent energy density with a typical 
scale which is large fraction of the typical eddy scale (Meneguzzi et al. 1981). Evidently, 
the random turbulent motions within eddies are capable of producing large fluctuations of 
the magnetic helicity on eddy scales. 

Assuming that this is the case, we can estimate the long wavelength tail of the 
magnetic field power spectrum by balancing the systematic dissipation of magnetic helicity 
inhomogeneities with their random generation via equation (11). (The dissipative term in 
this equation is not immediately apparent because it is divided between the two magnetic 
helicity current terms on the right hand side.) Given a typical eddy scale, A, with a velocity, 
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v, and fluctuating magnetic field, 6, we get a linear growth rate for the mean square 
magnetic helicity on some larger scale A as follows: 



where the factor of (A/A) 3 comes from considering the fraction of phase space containing 
the large scale modes. Since the dissipation rate is just A~ 2 (v 2 )t c , we find that 



In other words, on large scales the magnetic helicity will settle into a Poisson distribution 
after a dissipation time for the scale in question. This implies that 

/ \ \ 5 /4 



which is small, but not necessarily negligible. 

Let's consider instead the situation where the magnetic field is in equipartition with 
the turbulence on some very small scale, much smaller than the typical energy-bearing eddy 
size. That is, we consider A and A both small. In this case we need to replace the turbulent 
diffusion coefficient ~ (v 2 ) 1 / 2 ^ with ~ (v 2 ) l l 2 \A. Assuming a Kolmogorov spectrum for 
(v 2 ), this implies a turbulent diffusion coefficient which is larger by a factor (A/A) 4 / 3 . 
Nevertheless, equation (35) is still a fair estimate of the strength of the magnetic field on 
the scale A induced by random anomalous magnetic helicity currents, since we should also 
replace (v 2 )\ on the right hand side of equation (34) with (v 2 )\, which is larger by the same 
factor. In addition, we need to consider the action of scales intermediate between A and 
A, where the magnetic field power spectrum is described by equation (35). As we consider 
larger values of A the term (b 2 ) 2 in equation (34) declines as A~ 5 . On the other hand, since 
these scales are part of the inertial range of hydrodynamic turbulence, the factor r 3 scales 
as A 2 . Including the factor of (A/A) 3 as well we see that the magnetic helicity currents 




(33) 




(34) 




(35) 
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driven on all intermediate scales, and the self-interaction of fields on the scale A, are as 
important as the small scale driving force. We conclude that the field strength will increase 
exponentially at a rate which is comparable to the eddy turn over rate on the scale of 
equipartition, until the peak in the magnetic power spectrum approaches the peak in the 
turbulent kinetic energy power spectrum. This is consistent with the numerical simulations 
cited above, and lends support to the assumption of approximate equipartition between the 
magnetic field and the velocity field in a highly conducting turbulent medium. 

This has some interesting implications for the early evolution of the galactic magnetic 
field. If we imagine that the early galactic disk was turbulent in roughly the same way the 
it is today, that is with a typical fluid velocity of 10 km/sec, on scale ~ 100 pc, then we 
expect a disordered, equipartition magnetic field to be present after a few times 10 7 years. 
The diffusion time for length scales of a few hundred parsecs is an order of magnitude 
larger than the eddy turn over time at 100 pc, or about 3 x 10 8 years, at which point we 
expect the long wavelength tail to extend to magnetic domains with sizes comparable to a 
disk thickness. We are most concerned with annular domains, which resist shearing and 
can serve as the basis for generating a globally organized galactic magnetic field. These are 
oddly shaped, and we need to rewrite the scaling law in equation (35) as 



where N eddy is the number of independent turbulent eddies in a magnetic domain. This 
implies an annular magnetic field of about 10 _1 &, or ~ 3 x 10~ 7 gauss. The galactic dynamo 
operates at an efficiency which is some fraction of the rotation speed divided by the galactic 
radius, or a fraction of 10 _15 sec _1 . Even if the efficiency is as low as 10% the large scale 
field will reach equipartition in less than 10 9 years. The typical domain sizes should also 
continue to increase, but a discussion of that process is outside the scope of this paper. 

We stress that while the physical processes invoked in this argument are new, the basic 
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idea of a self-generated turbulent field as a critical step in the growth a galactic magnetic 
field is not new. The idea of generating a disordered field via turbulent fluctuations can 
be found in Kulsrud, Cen, Ostriker & Ryu (1997), although there the intent was to do 
without a large scale dynamo altogether. The physics of random fluid helicity fluctuations 
in a turbulent shearing medium has been previously discussed in Vishniac & Brandenburg 
(1997). In addition, the early growth of a galactic seed field, using the cross-helicity effect, 
has been discussed by Brandenburg & Urpin (1998) and Blackman (2000). This model 
differs in producing exponential growth of an arbitrarily weak small scale field using dynamo 
processes which respect magnetic helicity conservation. This growth proceeds through a 
kind of two-step inverse cascade. Large scale fluctuations in magnetic helicity are generated 
by the spontaneous appearance of regions of positive and negative magnetic helicity. These 
magnetic helicity densities are then transfered to large scale structures through dissipative 
coupling between the scales. 

4. Discussion and Conclusions 

We have examined the suppression of the dynamo effect due to the conservation of 
magnetic helicity. This has a substantial effect on the viability of dynamos operating in 
a closed box filled with homogeneous turbulence, the l a 2 dynamo'. Whether or not is 
completely suppresses them depends on whether or not there is a version of our model 
which generates a largely force-free dynamo field. On the other hand, it has remarkably 
little effect on the generation of large scale magnetic fields in differentially rotating systems. 
Dimensional estimates of dynamo growth rates are the same; the need for fully three 
dimensional turbulence to drive dynamo activity is unchanged; and the important role of 
differential rotation in producing large scale fields is unchanged. Quantitatively, we can see 
some differences. The most obvious one is that a previously ignored property of anisotropic 



-20- 



turbulence plays a critical role in dynamo activity and that there is a strong tendency 
for this quantity to have the correct sign for promoting dynamo activity. A more subtle 
effect is that the symmetry breaking requirements for a successful dynamo are significantly 
reduced. Finally, we note that a successful dynamo appears to rely on a large scale, spatially 
constant, magnetic helicity current. This current is not tied to the magnetic domain size, 
and its sign depends only on the sense of the differential rotation. 

We have also examined the role of random helicity currents in generating small scale 
magnetic fields. We find that this process plays a key role in allowing magnetic fields to 
reach equipartition levels in a turbulent medium. Using plausible numbers for turbulence 
in the galactic disk, we find that large scale magnetic fields can be generated at levels of 
~ 10 _1 of current values within the first few 10 8 years of the disk's existence, starting from 
very small scale, highly disordered and weak initial fields. The strength of the galactic scale 
seed field is irrelevant. 

It seems odd that eddy scale motions can generate significant large scale magnetic 
helicity currents when only a small amount of the large scale helicity is contained in the 
eddies. However, while the eddies cannot store any significant amount of magnetic helicity 
of one sign, they will create local fluctuations in the magnetic helicity in the presence of 
a large scale field, even when there is no large scale magnetic helicity. The typical size of 
these fluctuations will be ~ Ba ~ B 2 vt c . As long as these fluctuations can be moved in 
different directions, depending on their sign, there is no paradox inherent in our results. In 
particular we note that the anomalous magnetic helicity is of the same order as v times the 
typical amplitude of these fluctuations. 

We also note that while the usual mean-field dynamo violates magnetic helicity 
conservation for the large scale field, this is not equivalent to violating magnetic helicity 
conservation altogether. The total magnetic helicity that can be stored in individual eddies 
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is down from the magnetic helicity contained in a large scale field with the same amplitude 
by a factor of \ C /L. The implication is that the back-reaction which suppresses S m f ■ B 
won't set in until 



This should be unimportant when a fast dynamo of the kind discussed in the second section 
of this paper is operating. However, it may be important otherwise, both in astrophysical 
objects and in computer simulations. Of course, if a successful dynamo occurs, and creates 
a significant magnetic helicity density, then this may generate secondarily a kinetic helicity. 
We do not expect this effect to play an important role in the dynamo process. 

Finally, we have ignored the possibility that boundary conditions might play a role in 
the evolution of the magnetic field. Blackman and Field (2000a, b,c) have suggested that the 
ejection of magnetic flux across system boundaries could be a necessary part of the dynamo 
process (see also Kleeorin et al. 2000). This would appear to conflict with the results of 
'shearing-box' simulations of accretion disks (e.g. Hawley, Gammie & Balbus 1996, and 
Hawley 2000) which show the generation of a large scale field from initial conditions with 
no net magnetic flux and little large scale structure, even with periodic vertical boundary 
conditions. Our results confirm that the ejection of magnetic helicity is not a requirement 
for a successful dynamo. However, our results do imply that a large scale magnetic helicity 
current is a necessary part of the dynamo process, and in real systems (which typically lack 
periodic boundary conditions) this will usually lead to a magnetic helicity current across 
the system boundaries. A non-periodic computer simulation with closed boundaries will 
start to quench dynamo growth in a box of size L when the anomalous magnetic helicity 
current ~ B 2 (v 2 )r c , or equivalently ~ B 2 Dt, that accumulates in a magnetic domain at 
the boundary becomes comparable to the magnetic helicity such a domain can contain, or 
~ B 2 L B . This is a quenching rate of D T /L 2 B , that is, the turbulent diffusion rate for a 
magnetic domain. In a real system the magnetic helicity will be ejected in a wind. If the 
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wind dimension is Lw (presumably roughly the size of the system) then since the ejected 
magnetic helicity flux must be similar to the magnetic helicity flux required to run the 
dynamo we expect an energy flux in the wind from a slow rotator of at least 



where we have used equation (27). This energy flux is less than the energy available from a 
single magnetic domain in the dynamo system by the ratio of the eddy scale to the wind 
scale. The magnetic flux loss is bounded by B W V W) or 



For a slow rotator this limit can be more than the magnetic flux generation rate within a 
single magnetic domain. The implication is that either the wind contains flux elements of 
mixed sign, or the wind is slow compared to the escape velocity, or the wind becomes fast, 
but only at a distance comparable to V esc /Q. 

We are happy to acknowledge helpful discussions with Amitava Bhattacharjee, Eric 
Blackman, Patrick Diamond, Eun-Jin Kim, and Alex Lazarian. 



Rather than assume that all long wavelength variations in the magnetic helicity are 
equivalent, we could keep track of the magnetic helicity due to large scale field components, 
H = A ■ B, and that due to small scale field components, h = (a • b), separately. In this 
case, Jh is clearly associated with small scale structures, since it arises from eddy scale 
fluctuations in the electromotive force. We can write the equations for h and H as 
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A. The Transfer of Magnetic Helicity Between Scales 



d t h = -V • J H - 2B • (v x b) 



(Al) 
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and 

9 t iJ = V-((vxB)xA})+2B-(vxb) ( , (A2) 

where the subscript 't' denotes effects which transfer magnetic helicity between scales rather 
than between locations, and the sum of equations (Al) and (A2) is equation (11). The 
standard mean field electromotive force parallel to B falls into the category of scale transfer 
effects. However, it is calculated by neglecting any information flow from small scales to 
large, and consequently is not the entire term. Here we will show that the presence of a 
non-zero h leads to the efficient transfer of magnetic helicity between scales, which justifies 
lumping H and h together. 

If h is non-zero, then there is a piece of a which is correlated with b. We can write this 

as 

b 

a c = h—, (A3) 

where we adopt the same gauge choice (V • a = 0) used elsewhere in this paper. This in 
turn implies a correlation between the small scale magnetic field and the current, 

where we have used the fact that the length scale for h is much larger than an eddy scale. 
This extra piece of the magnetic field produces a first order change in the turbulent velocity 
field 

^ (B-V)Vxb 

Vc ~ 4np(V) kTc > (A5) 

which in turn produces a magnetic helicity transfer term of 

hr 

-2B-(v x b) t » -2— ^-B-((B • V)bx(V x b)>. (A6) 
Aiip{b z ) 

If the turbulence is approximately random and isotropic, then this reduces to 

-2B-(v x b) t » -2hr c V 2 A k\^ » -thr c k 2 Vl (A7) 
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where the directions _L and || are defined relative to the direction of the large scale field. 

Equation (A7) implies that small scale eddies dump their average magnetic helicity 
into large scale structures at a rate which approaches the eddy turn over rate as the large 
scale field approaches equipartition with the turbulence. This expression does not include 
the back reaction of the large scale field, which will be important if it already contains as 
much magnetic helicity as possible on energetic grounds. We will ignore this point here. It 
may have a quantitative effect on the saturation level of the magnetic field, although it is 
unlikely to be more important than turbulent dissipation. The transfer of magnetic helicity 
to larger scales can be treated as arbitrarily fast provided that the dynamo growth rate, 
F D , satisfies the inequality 

T D < r c (W A ) 2 ~ ^r-\ (A8) 

The implication is that for sufficiently weak large scale fields the dynamo will start off 
by creating large scale gradients in the small scale magnetic helicity, and subsequently 
evolve to the point where the magnetic helicity is transferred to larger scales as fast as it 
accumulates. This will set in when the transfer term given in equation (A7) is roughly 
equal to V • Jh or 

where L B is the length scale of the large scale magnetic field. This limit is the saturation 
value for the magnetic helicity on a scale k^ 1 times (/cL B ) _1 . In other words, the distinction 
between h and H will become dynamically unimportant long before h approaches its 
saturation value. 

In the presence of a non-zero fluid helicity we will also have the transfer of magnetic 
helicity between scales following the standard mean field dynamo formalism. This rate is 
comparable to V • Jh if we identify the symmetry breaking scale for the fluid helicity with 
the length scale of the large scale magnetic field. Consequently, when the large scale field 
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is weak enough that the transfer of small scale helicity to large scale helicities constitutes 
an impediment for the dynamo process, the standard mean field dynamo will be operating. 
Conversely, once the large scale field is strong enough to allow for efficient transfer of 
magnetic helicity between scales, the fluid helicity is a subdominant effect. 

Finally, in the absence of any large scale field at all, turbulent dissipation alone will 
transfer magnetic helicity to large scales by wiping out the smaller scale structures that 
generated the large scale variations in magnetic helicity. This process will move helicity 
to larger scales, at the dissipation time scale of the structures that contain the magnetic 
helicity fluctuations, until the containing structures are as large as the length scale for the 
magnetic helicity fluctuations. Since this will go to completion on the dissipation time 
scale for the largest scales involved, it is, by definition, too slow to dominate in successful 
dynamos. However, it will be important in the kind of stochastic processes described in 
section 3. It dominates over the helicity transfer rate given in equation (A7) when 




B 



(A10) 



or 




(AH) 
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